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ABSTRACT 

Aims. We study the distribution of the photometric rotation period (Prot), which is a direct measurement of the surface rotation at 
active latitudes, for three subsamples of Sun-like stars: one from CoRoT data and two from Kepler data. For this purpose, we identify 
the main populations of these samples and interpret their main biases specifically for a comparison with the solar P^oi- 
Methods. Py^t and variability amplitude (A) measurements were obtained from public CoRoT and Kepler catalogs, which were 
combined with public data of physical parameters. Because these samples are subject to selection effects, we computed synthetic 
samples with simulated biases to compare with observations, particularly around the location of the Sun in the Hertzsprung-Russel 
(HR) diagram. Publicly available theoretical grids and empirical relations were used to combine physical parameters with Piot and A. 
Biases were simulated by performing cutoffs on the physical and rotational parameters in the same way as in each observed sample. 
A crucial cutoff is related with the detectability of the rotational modulation, which strongly depends on A. 

Results. The synthetic samples explain the observed distributions of Sun-like stars as having two main populations: one of young 
objects (group I, with ages younger than ~ 1 Gyr) and another of main-sequence and evolved stars (group II, with ages older than ~ 
1 Gyr). The proportions of groups I and II in relation to the total number of stars range within 64-84% and 16-36%, respectively. 
Hence, young objects abound in the distributions, producing the effect of observing a high number of short periods around the location 
of the Sun in the HR diagram. Differences in the Pyot distributions between the CoRoT and Kepler Sun-like samples may be associated 
with different Galactic populations. Overall, the synthetic distribution around the solar period agrees with observations, which suggests 
that the solar rotation is normal with respect to Sun-like stars within the accuracy of current data. 

Key words. Stars: rotation — Stars: evolution — Stars: Solar-type — Sun: rotation 


1. Introduction 


A question of high interest is how exactly the Sun will evolve 
and, as a consequence, how it will affect the planetary envi¬ 
ronment. Although the Sun’s evolution has long been modeled 
from thermonuclear reaction rates, explaining non-canonical ef¬ 
fects, such as sunspots and magnetic fields, is still puzzling. 
Clearly understanding the solar structure and evolution is nec¬ 
essary to know, for example, how long our central star will 
maintain the conditions required for life in the solar system, 
or whether its magnetic activity can damage our technologi¬ 
cal needs, among other concerns (e.g., Melott^Diom^ 201 ll: 
iJonesI I 2 OO IShibata et al.ll2013l ISteinhilber & Beej 120131) . To 
accurately predict the near and distant future characteristics of 
the Sun, it is essential to understa nd its past and project its 
future by studying other stars (e.g.. ICharbonneaull2010[ 120131: 


* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with the participation of ESA’s RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. 


lEkstrom et al.ll201^ . Specifically, stars with a similar mass and 
chemistry to that of the Sun may provide important constraints 
for th e theoretical models used in predicting the Sun’s history 
(e.g., ISchuler et akl 1201 It iDatson et akl l2012l) . To refine these 
models, stellar rotation is fundamental because this parameter is 
a key feature in controlling the root-cause of the s tructure, evolu¬ 
tion, chemistry, and mag netism of the stars (e.g.. lPalaciosll201^ 
iMaeder & Mevnenl2014l) . In particular, the measure of the sur¬ 
face rotation can substantially improve observational constraints 
on theory from a proper modeling of the stellar i nterior (e.g., 
IZahnlflQ^ lEkstrbm et al.ll20f3 iBrun et aTll2014l) . Eurther re¬ 
finements can be obtained from the internal rotation, which has 
recently been measured in the Sun and several different types of 
stars thanks to remarkable advances in a steroseismology (e.g., 
iGoupil et al.ll2013T:lDeheuvels et al.ll2014l) . 

Another key question is how similar or different the Sun is 
compared to Sun-like stars, and rotation may also be crucial i n 
this comparison (e.g.. lGustafssonl[r998l: lGonzalezlll99^ 1200ll) . 
Different studies have discussed whether the Sun rotates nor- 
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mally with respect to its Sun-like count erparts, at least co nsid¬ 
ering the su rface rotation. For example, ISoderbloml (Il983h and 
iGravl (Il982h suggested th at the Sun rotates no rmally for its age, 
a finding corroborated by iRobles et ^ (l2008h . However, as re¬ 
ported by Ide Freitas et al.l (12013 ), these studies were based on 
the projected rotational velocity (v sin /), which is an indirect ro¬ 
tation measurement wit h the intrinsical l y unkn own axis orienta¬ 
tion angle i. In contrast, iMetcalfe et al.l (12013h found in a recent 
study of magnetic activity cycles of the solar analog e Eri that 
this star undergoes long-term magnetic cycles compatible with 
those of the Sun, but the rotatio n period Prnt deter mined from 
the photometric time-series ('e.g., IStrassmeieii[2009l) is approxi¬ 
mately twice as fast. Based on Pi-ot, which is a direct measure¬ 
ment of the surface rotation at the active latitudes, it is suggested 
that the Sun may rotate more slowly than its analogs. Neverthe¬ 
less, the question about solar rotation normality remains unan¬ 
swered, especially when we consider its rotation compared to 
direct stellar rotation measurements. 

To date, most comparative studies of the Sun’s rota¬ 
tion relative to other stars ar e based on the projected rota¬ 
tional velocity from v sin i (e.g.- lPorto de Mello & da Silvall99'^ 
iMelendez et al.ll^06h , which, as explained, may hinder proper 
comparison. In contrast, the number of databases with direct 
Frot measurements from the photometric ti me-series is grow¬ 
ing rapidly (e.g., Strassmgieret ^ l200d: IStrassm eied 120091: 
Hartman_et^ 20 id: De Medeiros et al.l 120131 : iReinhold et^ 


2013t McQuillan et al.l I 2 OI 4 I 1 . Current examples of databases 


with stellar rotation periods are th ose obtained from the CoRolQ 
(iBaglin et alj|2009h and Keple^ (iKoch et alJl2010l) space tele¬ 
scope data, which include long-term light curves (LCs) with 
high-temporal resolution and high-photometric sensibility for 
more than 300,000 objects. Together, these databases include 


or candidates (e.g.. 

Meibom et alj|201 It 

Affer et al.ll2012l 2013; 

De Medeiros et alJ 

2013; Nielsen et al.l 

2013t McOuillan et alJ 

201 3allhll2014HReinhold et alJl201 3t Walkowicz, & Rasrill201 3l). 


However, for many of these stars, no spectroscopic data are 
available. Thus, it is difficult to associate their periods with other 
physical parameters, such as metallicity, mass, radius, tempera¬ 
ture, and age. 

In particular. lMcOuillan et al.l (I2013ail2014ll found evidence 
of a bimodal Frot distribution for cool stars that becomes shal- 
lo wer with higher t emper atures. Such a bimodality, also detected 
in IReinhold et al.l (12013h . has been proposed to originate from 
stellar populations that evolved differently from one another. The 
current lack of spectroscopic data is one of the main limitations 
in explaning this effect. The Kepler database so far provides pho¬ 
tometric estimates of physical parameters for nearly all targets 
that are useful for several purposes. With regard to CoRoT, these 
photometric estimates are typically available for the targets of 
the asteroseismology channel. Thus, follow-up is crucial to ob¬ 
tain this physical information for the targets of the exoplanet 
channel (i.e., those considered in the present work). 

Recently, iGarcfa et al.l (12014ll examined the rotational evo¬ 
lution of solar-like stars based on the Kepler data. These au¬ 
thors performed a refined study of the age-rotation-activity re¬ 
lations by combining rotation period measurements with pre¬ 
cise asteroseismic ages. One of the results was an analysis of 
the evolution of the rotation period as a function of stellar 
age for selected 12 cool dwarfs. The period-age relation corre¬ 
sponded well with previous calibrations, such as the Skumanich 


' http://smsc.cnes.fr/C0R0T/ 
^ http://kepler.nasa.gov/ 


law iSkuriiunich[ 19721) and t hose o btained by iBarne i (120071) 
and iMamaiek & Hillenbrandl (l2008l) . Based on the referenced 
work, a fit with the 12 selected stars is compatible with the Sun 
location in the period-age diagram. However, the Sun lies at a 
slightly longer period than the fit, a small discrepancy that needs 
to be investigated in more detail. 

For general studies of Kepler stars, iHuber et"^ (1201 4l) pro¬ 
vided a compilation of physical parameters obtained from differ¬ 
ent methods with improved values and well-computed errors by 


adjusti ng observational da ta with theoretical grids. For CoRoT 
targets, ISarro et ^ (1201 3l) recently published a relatively long 
list of physical parameters obtained automatically from spec¬ 
troscopy. These data represent a unique opportunity for follow¬ 
up of the photometric periods determined from CoRoT and Ke¬ 
pler observations. In addition, a comparison between observa¬ 


tions and theory can be performed by considerin g recent evolu¬ 
tionar y models developed by the Geneva Team (lEkstrom et al.l 
I 2 OI 2 I) . which account for stellar rotation based on a rich set of 
physical ingredients (see Sect. O. Accordingly, photometric pe¬ 
riods can be compared with these models. It is necessary for such 
a comparison to carefully analyze the observational bias of the 
observed sample, as described below. 

Biases are typically modeled in st ellar population sy nthe- 
sis, such as in the TRILEGA^ (e.g.. iGirardi et al.]l2005l) and 
SYCLISIQ (iGeorgy et al.ll201^ codes. TRILEGAL combines 
nonrotating evolutionary tracks with several ingredients, such 
as the star formation rate, age-metallicity relation, initial mass 
function, and geometry of Galaxy components, to simulate, in 
particular, the distribution of physical parameters and photomet¬ 
ric measurements of different stellar populations. SYCLIST is a 
recent code developed by the Geneva Team that includes stel¬ 
lar rotation in a detailed set of ingredients to compute synthetic 
parent samples. The code considers anisotropies of the stellar 
surface, such as the latitude-dependence of temperature and lu¬ 
minosity and their effects on the limb darkening, to predict how 
the measured physical parameters can be affected by different 
rotation axis orientations. The resulting biases can, for example, 
affect the age estimation of coeval samples or of individual field 
stars because the physical parameters may deviate from their ac¬ 
tual values. Although the SYCLIST code can be used to predict 
deviations on rotation measurements that are the result of differ¬ 
ent axis orientations, it does not yet consider the detectability of 
Fiot measurements from EC variations. 

Eor a deeper study of stellar rotation, understanding biases 
of Frot measurements can be particularly useful. As such, an 
important effect must be considered: the amplitude of the ro- 


tational modulation decreases as long as the stars evolve (e.g ., 
iMessina et alJl200il 120031 : lGiampapal201 lb ISpada et ^120111) : 


this decreases their detectability. As a consequence, current field 
stellar samples with Fjot measurements are biased, exhibiting a 
lack of stars around the solar rotation periocQ. This bias hampers 
a proper analysis of the solar rotation normality regarding the 
Frot distribution of Sun-like field stars. Therefore, a quantified 
analysis is needed to better identify these biases in field stellar 
samples with variability measurements. 

We here study the empirical distribution of Frot for field stars 
with physical parameters similar to those of the Sun. To this 
end, three different subsamples of Sun-like stars were obtained 
from the CoRoT and Kepler public catalogs, and synthetic sam- 


^ http://stev.oapd.inaf.it/cgi-bin/trilegal 
http: //ohswww. unige. ch/Recherche/evoldb/index/ 

^ The solar period ranges from 24.47 days at the equator to 33.5 days 
at the poles and has an avearage of 26.09 days dLanza et alJl200.3l) . 
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pies were built by reproducing the main biases of the observed 
samples. Actual and synthetic samples can then be compared to 
one another to interpret the observational biases and to study 
the solar rotation normality. This analysis is performed by test¬ 
ing the most recent theoretical predictions combined with a large 
amount of new observational data. Combining this observational 
and theoretical information in the same analysis allows us to 
conduct an unprecedented study concerning these CoRoT and 
Kepler targets compared with the solar rotation. 

We begin by describing our methods (Sect. ^ for sample 
selection from the public CoRoT and Kepler data. Next, we 
explain the methods we used to build the synthetic samples 
(Sect. [3). We discuss the rotation period distribution of young 
stars (Sect. S, which were needed to constrain our synthetic 
samples. We present our results (Sect. IS starting from the val¬ 
idation between the observational data and theoretical predic¬ 
tions for the Sun-like CoRoT sample. This sample, obtained 
from iDe Medeiros et alJ (1201 3l hereafter DM13) was carefully 
considered as a catalog of rotating candidates because physical 
information of its sources was very limited during its develop- 
ment. Now , new s pectroscopic physical parameters provided by 
ISarro et ^ (1201 3h allow us to determine whether this sample is 
indeed a catalog of Prot- Finally, the Prot distribution is analyzed 
for the CoRoT and Kepler subsamples of Sun-like stars by inter¬ 
preting biases and by comparing the stars with the solar rotation. 


2. Sample selection and physical parameters 

In this section, we summarize the catalogs used to obtain our ob¬ 
servational samples, from which we constructed the correspond¬ 
ing synthetic samples. These include the CoRoT and Kepler cat¬ 
alogs with variability period measurements and those used to 
derive the physical parameters. 


2.1. Selection of CoRoT targets and their parameters 


ulation caused by dynamic star spots (semi-sinusoidal variatiorQ 
def ined in DM 1 3), and the final parameters were refined. 

ISarro et ^ (1201 3li used neural networks similar to those 
employed in the CoRoT Variable Classifi er (C VC) tool 
(iDebosscher et al.l l2007l 120091: iBlomme et alJ l20ld) to obtain 
physical parameters from spectroscopic data of 6,832 CoRoT 
targets. The method was applied to FLAMES/GIRAFFifl spec¬ 
tra by using two different types of training sets. One is com¬ 
posed of synthetic spectra from Kurucz models (iBertone et al.l 
1200^ and TLUSTY grids (iHubenv & Lanzl[l9^ (KT set) and 
the other of actual ELODm spectra with known parameters 
from the ELODIE library ( Prugniel & Soubir^ 12001 1 l2004h 
(ELODIE set). Physical parameters, particularly gravity (logg) 
and effective temperature (Teff), were then obtained indepen¬ 
dently for both sets. The main limitations of this catalog are the 
uncertainties in the computed parameters. Thus, no metallicity 
values are available and the errors in logg and Teff are relatively 
high (originally with typical values of approximately 0.5 dex and 
400 K, respectively). As a first step to optimize our analysis de¬ 
spite the limitations, we averaged multiple independen t estima¬ 
tions of log g and T^g provided by ISarro et ^ (1201 3h for each 
target, which decreased the typical uncertainties slightly to ap¬ 
proximately 0.4 dex and 300 K, respectively. 

Eor the sample selection, we combined the catal og of DM13 
with t he available physical parameters given by ISarro et ^ 
(l2013h . This provided a subsample of 671 targ ets, namely 
the DM S sample, which we analyze here. Because ISarro et 
(1201 3l) did not provide metallicity values, we assumed that the 
sample has a mixture of metallicities with an average around the 
solar value. This assumption is reasonable given that the over¬ 
all metallici ty \MfH] distr i bution in the CoRoT fields, in accor¬ 
dance with iGazzano et al.l (1201 Oli . has a mean value and stan¬ 
dard deviation of -0.05 + 0.35 dex. Eor a particular analysis of 
stars with T^g and logg similar to the solar values, we selected a 
small subsample of 175 stars, namely, the DMS© sample, within 
an elliptical region of the size of typical uncertainties around 
Teg - 5772 + 300 K and logg = 4.44+ 0.4 dex. Then, we studied 
the distribution of Prot for this sample. Despite the large errors 
in log g and Teg and the assumptions on metallicity, the DMS© 
sample can be interpreted from a comparison with a simulated 
sample, which considers the main observational biases and un¬ 
certainties, as described in Sect. [2 


Erom a parent sample of 124,471 LCs, DM13 presented a cat¬ 
alog of 4,206 CoRoT stellar rotating candidates with unam¬ 
biguous variability periods. A detailed procedure was developed 
for proper selection of the targets, based only on CoRoT LCs 
and 2MASS photometry, before additional data such as spec¬ 
troscopy data were available. The procedure includes an auto¬ 
matic pre-selection of LCs that corrects them for jumps, long¬ 
term trends, and outliers, and identifies a preliminary list of peri¬ 
ods, amplitudes (A), and signal-to-noise ratios (S / N). The vari¬ 
ability period was computed from Lom b-Scargle (lLomblll976l : 
IScargleiri982t [Home & Baliunasl[l98^ periodograms, the am¬ 
plitude from a harmonic fit of the phase diagram for the main 
period, and the S /N was the amplitude-to-noise ratio calculated 
from the EC. These parameters were then used to select the high¬ 
est quality LCs (high S /N) in the expected amplitude and period 
ranges for rotating candidates. This criterion is related to the de¬ 
tectability of rotational modulation in an EC and was used in 
the present work as a relevant ingredient to produce a synthetic 
sample with biases (see Sect. B- Einally, visual inspection was 
performed by applying a list of well-defined criteria for identify¬ 
ing photometric variability compatible with the rotational mod- 


2.2. Selection o/Kepler targets and their parameters 


Several recent works (iM cOuilla n et al.l I20l4 l2013allbt 
Nielsen et al.l 120131: IWalkowicz & Bas^ 20131: iReinhold et al.l 
20131) provide rotation periods for a lyge por tion of the Kepler 


20131) provide rotation periods for a lyge por tion of the Kepler 
targets. In particular, IReinhold et al.l (l2013l) derived periods 


for 24,124 targets observed during quarter Q3 by following an 
automatic method with partial inspection. In the method, the 
range between the 5th and 95th percentile of the nor malized flux 
distribution of a EC, named variability range Rvar (iBasri et al.l 
l2()Io[l20ni) . is an adimensional quantity with a numerical value 


^ Semi-sinusoidal variability is characterized by some asymmetry of 
the maximum and minimum flux with respect to its average over time, 
somewhat irregular long-term amplitude variations, and semi-regular 
multi-sinusoidal short-term flux variations, typically with an amplitude 
< 0.5 mag and period > 0.3 days. This description is based on the 
dynamic behavior of star spots, as observed for the Sun (see DM13, 
Sect. 2.2). 

’ http://WWW.eso.org/sci/facilities/paranal 
/instruments/flames/ 

* http: //WWW . ohs-hp. fr/guide/elodie/elodie-eng . html 
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approximately equal to the variability amplitude A in units of 
magnitude. The detectability of rotating stars was primarily de¬ 
termined by defining an f?var threshold (of 3%c), which produced 
a bias concept ually comparable to the S/N threshold of DM13. 
The catalog of iReinhold et al.l (12013h provides rotation periods 
ranging from 0.5 to 45 days, thus covering a region that can be 
compared with the solar rotation pe r iod. 

More recently. iMcOuillan et al.l (l2014l) provided a new cata¬ 
log of rotation periods ranging from 0.2 to 70 days for 34,040 
Kepler targets observed during quarters Q3-Q14 (which is a 
much longer time span than considered in IReinhold et al.ll2013h . 
This is currently the largest Kepler sample of rotation periods 
and was obtained from a fully automatic method named Au- 
toACF. The method is based on the autocorrelation function, 
which was applied to a trained neural network that selected a 
sample of stars exhibiting rotational modulation. For the pe¬ 
riod detectability, the authors defined a weight w that has a so¬ 
phisticated calculation based on the ACF periodogram, physical 
parameters, and neural network selection. This weight, with a 
threshold Wthres = 0.25, produces a bias in the detectability of the 
rotational modulation, somewhat similar to the aforementioned 
S I N or Ryar parameter s. 

iHuber et’aH (1201 4l) recently presented a catalog of physical 
parameters for Kepler stars that can be combined with rotation 
period measurements. These authors provided a large compila¬ 
tion of physical parameters estimated from different methods 
for 196,468 stars observed in the Kepler quarters 1 to 16. For 
the majority of the targets, atmospheric properties (temperature, 
gravity, and metallicity) were estimated from photometry and 
have high uncertainties. In addition, there is also a noticeable 
number of objects with more refined measurements obtained 
from spectroscopy, asteroseismology, and exoplanet transits. Fi¬ 
nally, the authors refined the parameters by adjusting observa¬ 
tions with theoretical grids. For this reason, this catalog has a 
sharp distribution in the Hertzsprung-Russel (HR) diagram, as 
shown in Sect. |2 which is also suitable for our methods used 
to produce a synthetic sample. On average, the typical uncer¬ 
tainties in Teff, logg, and metallicity {FejH} lie in the refereed 
catalog at approximately 170 K, 0.2 dex, and 0.2 dex, respec¬ 
tively. To obtain a subsample of Kepler Sun-like stars, we se¬ 
lected only targets with metallicity -0.2 dex < {FejH} < 0.2 dex 
and defined an elliptical region in the HR diagram within - 
5772 + 170 K and logg = 4.44 + 0.2 dex . This region w a s used 
to s elect Sun-like stars in th e catalogs of IReinhold et al. (201_3h 
and McQuillan et ^ (l2014t) . each combined with iHuber et all 
(l2014li . namely, the RH^t, and MHq samples, which are composed 
of 1,836 and 2,525 objects, respectively. The same metallicity 
range and HR diagram region was used to produce the synthetic 
forms of these two Kepler Sun-like samples. 

3. Synthetic sample and main biases 

lEkstrom et al.l (1201 2l) presented a new version of the Geneva stel¬ 
lar evolution code dEggenberger et al.ll200^ for noninteracting 
stars that includes theoretical predictions for the surface rotation. 
To predict this parameter, the authors considered several physical 
assumptions regarding the stellar interior. T he main basi s of this 
model is the shellular-rotation hypothesis (IZahnlflOQ^ . which 
states that the turbulence is substantially stronger in the horizon¬ 
tal than in the vertical direction, yielding nearly constant angular 
rotation in each isobaric region. On this basis, horizontal dif¬ 
fusion is combined with meridional circulation and shear turbu¬ 
lence to predict transport mechanisms of angular momentum and 
of chemical species. The angular momentum is conserved during 


the whole stellar evolution, considering the mass-loss influence. 
Einally, both atomic diffusion and magnetic braking in low-mass 
stars are accounted for in a homogeneous fashion. The model 
provided grids of evolutions of solar-metallicity stars for masses 
ranging from 0.8 to 120 Mq, with low-mass stars evolving from 
the zero-age main sequence (ZAMS) to the core helium-burning 
phase. Therefore, these grids provide a theoretical distribution of 
Prot without observational biases. Erom this theoretical distribu¬ 
tion, synthetic biases and uncertainties can be applied to study 
how these effects may modify actual observations. Our analy¬ 
sis aims to resolve the actual distributions of three Sun-like 
samples (the DMS©, RHq, and MHq sample) b y an inversion 
proced ure, where we used the theoretical grids of Ekstrom et al.l 
(12012l) as a root ingredient to compute the three respective syn¬ 
thetic samples. Eor these calculations, we identified the main 
biases and defined appropriate parameters to describe each P^ot 
distribution. These parameters were then properly adjusted to fit 
the synthetic distributions with the observations. 

In summary, a synthetic parent sample of solar metallicity, 
obtained from an actual parent sample (see details below), was 
generated with the parameters Teff and log g assumed to be error- 
free. These data were used as input to add the following parame¬ 
ters to the parent sample: P^oi, A, S/N, and w, also assumed to be 
error-free. These parameters were set bas ed on a combination o f 
ingredients - such as theoretical grids of lEkstro m et al. (l2012h. 
expect ed P^ot evolution of young stars from iGal let & Bouvier 
(2013) , and P^t versus A empirical relations from iMessina et al. 
(120011120031) - as detailed below. Next, random fluctuations were 


applied to all the input and output parameters according to their 
errors in the actual sample. Erom this step, only virtual targets 
with a certain detectability parameter greater than a given thresh¬ 
old were selected. This parameter was defined according to each 
Sun-like sample to reproduce potential biases originating from 
the instrumental capability of detecting rotational modulation. 
Eor the DMSq sample, the S /N was used to perform a cutoff sim¬ 
ilar to that stated in DM13. Eor the RHct, sample, A was used t o 
mimic the selection by Rvar performed in IReinhold et al.l ( l2013l) . 
Eor the MHq sample, an empirical relation between A and w 
was identified and used to generate random synthetic values of w 
from the previously drawn A values. Eor every synthetic sample, 
we tested the bias produced by a cutoff with different parameters, 
S IN, A, and w, to discuss the effects on the Prot distribution (see 
Sect. 13. Einally, only stars with properties similar to the Sun 
were selected by applying the same boundaries as for each of 
the three actual samples. Because our final samples have phys¬ 
ical parameters similar to the solar values, the synthetic parent 
samples did not need to simulate stars well that are very differ¬ 
ent from the Sun. 

To further elaborate, the synthetic sample was based on the 
following assumptions: 

i. The parent distr ibution of each Su n-like sample is similar 
to the sample of Smo et akl (1201 3l) for CoRoT targets or to 
that of iHuber et al.l (1201 4l) for Kepler targets. 

ii. The rotation period follows a spread distribution for stars 
younger than ~1 Gyr and conv erges to a common evolu tion 
per stellar mass, as proposed bv iGallet & BouviCT(l2013l see 
their Eig. 3). 

iii. Based on the latter assumption, the synthetic parent sample 
is divided into two groups: one of young stars (group I), with 
ages younger than ~1 Gyr, and another of MS and evolved 
stars (group II), with ages greater than ~1 Gyr. Eor group I, 
the actual period distribution is difficult to determine accu¬ 
rately (see Sect. Si; thus, random values can be set within 
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a modeled distribution to analyze how they may be super¬ 
imposed with group II. For group II, the period follows the 
theoretical predictions provided bv lEkstrom et al.l (l2012l) . 
iv. The detectability of rotational modulation depends on a cer¬ 
tain parameter threshold: S IN for the DMSq sample, A for 
the RHq sample, and w for the MHq sample. 

V. The detectability parameter is strongly correlated with the 
variability amplitude or is the amplitude itsel f, which is re¬ 
lated to the rotation period, as reported by iMessina et al.l 

(l2nniL[2nn^ . 

vi. The final Prot distribution is also affected by uncertainties 
in Teff and log g, which cause the sample to become more 
widely distributed. 


To reproduce the parent distribution (as listed above in 
item i), w e first computed a H ess-like HR diagram of the original 
sample of jSarro^'^ (1201 3h for the CoRoT targets or a subsam¬ 
ple of iHuber et al.l (1201 4ll with metallicity -0.2 dex < [Fe/H] < 
0.2 dex for the Kepler targets. These distributions were then used 
as a map of probabilities, as shown in Fig.[T] The map resolution 
of the Kepler sample is higher than that of the CoRoT sample 
because of the considerably greater number of observed Kepler 
targets compared with CoRoT targets. Therefore, a more refined 
study can be developed from the Kepler samples. Actual samples 
were used as parent samples instead of population syntheses be¬ 
cause our methods aim to reproduce the detectability bias of the 
rotational modulation and do not account for how the parent pop¬ 
ulation was produced. As such, all biases related to the observed 
fields are already implicit in the considered parent samples. Pop¬ 
ulation syntheses can also be used in this method if desired. As 
another possibility, an image deconvolution could be applied to 
the parent sample probability map to reduce the error spread. We 
tested this potential effect us ing a Ri chardson-Lucy a lgorithm 
(iRichardsonI 1 972t iLucvl 1 97^ see also lLeao et al.l2()0^ and ob¬ 
tained similar results to those shown in Sect. |5] with some fluctu¬ 
ations. This suggests a robustness of the method, which produces 
stable results for different modelings of the parent distribution. 

From each of the maps depicted in Fig. [T] a synthetic parent 
sample with 10^ stars was generated randomly with probabilities 
respecting the map distribution. The drawn values of and 
logg were then defined as being exact for the synthetic sample. 
T hen, the genera t ed sam ple was plotted with the theoretical grid 
of lEkstrom et alJ (1201 2l) to perform the following definitions: 


i. Eor Teff out from the grid domain, the nearest valid value 
was set for a synthetic star. 

ii. Synthetic stars gene rated below the 1 Gyr isochrone of 
lEkstrom et alJ (l2012l) were defined to belong to group I, the 
remaining stars were designated to compose group II. 


This means that all stars generated from the parent sample prob¬ 
ability maps were used, including those drawn out from the grid 
domain. 

When groups I and II were established, theoretical predic¬ 
tions and empirical relations were used to aggregate Prot, A, and 
S /N values to the synthetic parent sample. Eirst, random periods 
were set to the synthetic stars of group I by following a custom 
distribution. The empirical distribution of this region is discussed 
in Sect. 0] to propose a synthetic distribution that may be super¬ 
imposed with group II. Eor the synthetic stars of group II, their 
Teft and log g values (initially defined as being exact) were used 
to set them to the co rresponding rotation periods given in the 
lEkstrom et alJ (1201 2l) theoretical grid. The grid, valid for the so¬ 
lar metallicity, was built by interpolating theoretical tracks and 
isochrones altogether. 
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Fig. 1. Probability distribution in the HR diagram for the parent sam¬ 
ples considered in this work. Top panel: distribution for th e CoRoT 
field, obtained from the full sample of ISarro et ^ ( l2013h . Bottom 
p anel: distribu t ion fo r the Kepler field, computed from a subselection 
of iHuber et all ( 120141) with a metallicity -0.2 dex < [FejH} < 0.2 dex. 
The red-green-blue color gradient in the map represents decreasing 
probability; the high end is th e red part of the gradi ent. Solid lines in¬ 
dicate the 1 Gyr isochrone of lEkstrom et all ( l2012h . and dashed lines 
illustrate the 1 Mq theoretical track provided by the same authors. The 
Sun is represented by its standard symbol. The ellipses around the so¬ 
lar location depict the subsamples of Sun-like stars (Sun-like samples). 
These were selected equally within the typical uncertainties in Teff and 
log g for each parent sample for observed and synthetic data (see text). 


To set up A values for the synthetic sample, we considered 
the empirical relations of Emt versus A obtained bv iMessina et al.l 
(I 2 OOII l2()()^ . which were sorted according to spectral type. In 
these relations, the stars are somewhat uniformly distributed be¬ 
low their corresponding A(Prot) empirical functions. The actual 
unbiased distribution in the P^ot versus A diagram cannot be de¬ 
rived from observations because observed data are biased. Such 
a distribution could be obtained from theoretical models if spot- 
induced variability amplitude calculations were included in the 
physical ingredients; however, such calculations were not imple¬ 
mented in the current models. Eor the synthetic stars, A values 
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were therefore set u p as uniform random value s below A(f’iot) 
curves (provided in iMessina et al.l [200111200^ . These curves 
were used for different spectral types, whic h were defined ac¬ 
cordi ng to the synthetic values of Teff (see iHabets & Heintzel 
Il981h : M-type for T^g < 3700 K, K-type for 3700 K < T^g < 
5200 K, and G-type for T^g > 5200 K. 

The approximation of considering a un iform probabil i ty dis¬ 
tribut ion below the empirical relations of IMessina et al.l (1200 ll 
l2003l) may affect the final period distribution, especially at short 
periods, where the period-amplitude diagram is more widely dis¬ 
tributed. One main effect in the simulations is a change in the 
proportions of groups I and II with respect to the total num¬ 
ber of stars. Pi and pu, respectively. These proportions are also 
sensitive to the line that defines these groups (here the 1 Gyr 
isochrone), which may fluctuate somewhat around its assumed 
HR diagram region and may have systematics. To adjust these 
limitations in the present work, the ratio pilpn was therefore 
set as a free parameter. Adjusting these proportions might be 
thought to produce a degeneracy: the synthetic Prot distribution 
would fit the observations without the need of a detectability 
threshold such as S fN, A, or w. We checked that combining this 
adjustment with a proper detectability threshold is certainly im¬ 
portant to obtain an optimal fit, otherwise, a noticeable discrep¬ 
ancy can occur, as explained in Sect. 15.31 For the CoRoT sample, 
a detectability threshold is mandatory to obtain an acceptable fit 
while it has a weaker effect on the Kepler samples, but it does 
contribute to improve the fits. Still, the Prot distribution is bi- 
modal regardless of its detectability bias. 

Furthermore, empirical magnitude-noise diagrams were built 
from actual data: one from CoRoT data and another from Ke¬ 
pler data (see Fig.|2]i. These maps were used to generate random 
magnitudes and noise levels for each corresponding synthetic 
sample. We estimated the noise level based on Eq. (1) given 
by DM13. This is a simple filtering of the high-frequency con¬ 
tribution that approximates the actual noise when the variability 
signal has a period considerably longer than the LC cadence. 
The noise level estimated from this equation has the advantage 
of being simplistic because no prior fit is needed for the calcula¬ 
tion. We verified that this noise level exhibits a similar behavior 
as a noise level estimated from the residual of a harmonic fit or 
a boxcar smooth of the LC. The two maps shown in Fig. |2]have 
two branches each, which are explained by the different types 
of LC bins (astero- versus exo-fields for CoRoT or short- ver¬ 
sus long-cadence for Kepler). Once noise levels were set to all 
synthetic stars, the amplitudes drawn above were used to set rea¬ 
sonable S IN values for the synthetic sample. In principle, S /N 
is not needed to simulate the Kepler samples, but it was used to 
test how the Pjot distribution depends on the cutoffs of different 
detectability parameters (see Sect. 1531) . 

Generating w values, which are needed to simulate the MHq 
sample, is particularly challenging because it is sophisticated. 
Lor an approximate solution, we considered an empirical dis¬ 
tribution of A versus w, which was used as a probabil ity map 
(see Lig.O. In the public data of lMcOuillan et al] (l20I4l) . w val¬ 
ues are available for the entire parent sample, but amplitudes 
are only given for w > 0.25. Hence, we computed indepen¬ 
dent A values from all Kepler LCs of the parent sample. This 
w as performed by automa tically combining multiple quarters, as 
in lMcOuillan et all (l2014ll . detecting the main periods of Lomb- 
Scargle periodograms and computing the amplitude of their cor¬ 
responding phase diagrams, regardless of their variability na¬ 
tures. This simple procedure produced amplitu de estimates that 
were h ighly compatible with those provided in lMcOuillan et al] 
(I 2 OI 4 I 1 . Lrom the probability map of Lig. |3l each synthetic w 
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Fig. 2. Probability maps for the empirical magnitude-noise diagrams of 
the CoRoT (top panel) and Kepler (bottom panel) parent samples. 


value was obtained by respecting the probability distribution of 
a section of the map according to a synthetic A value drawn pre¬ 
viously. 


Next, random fluctuations were applied to the synthetic 
generated values of T^^g, logg, Prot, A, and S jN according 
to their observational errors (no error was established for 
w). Then, a synthetic subsample with SIN, A, or w values 
greater t han a threshold wa s selec t ed as being a biased sam- 


ple, as in De Medeiros et al.l (l2013l) . iReinhold et al.l (12013l) . and 


iMcOuillan et al.l(12014 ). These thresholds were computed as free 
parameters to obtain the best fits between the synthetic and the 
observed Piot distributions. Linally, an elliptical region around 
the solar T^g and log g values was defined according to their typi¬ 
cal uncertainties, and a synthetic subsample of Sun-like stars was 
obtained. A particular cutoff of Prot < 45 d was also performed 
for the RH 0 sample and another of T^g < 6500 K was applied for 
the MHq sample to mimic similar cutoffs in the observed data, as 
performed in the original works. Best fits between synthetic and 
actual samples were computed from the Kolmogorov-Smirnov 
(K-S) test by minimizing the distance between their cumulative 
distributions. Uncertainties were estimated by analyzing the pa- 
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Variability amplitude (mmag) 


Fig. 3. Probability map for the empirical amplitude-weight diagram, A 
versu s w, of the Kepler paren t sample. Public w values were obtained 
from iMcOuillan et al.l ( |2014) . but we independently computed the A 
values. 


rameter ranges in which the K-S distance was smaller than twice 
its minimum. 


4. Period distribution of young stars 

iGallet & Bouvied (1201 3h developed models to predict the rota¬ 
tional evolution for pre-main sequence (PMS) stars. In these 
models, stars with different initial angular velocity values evolve 
by keeping this parameter constant during the disk accretion 
phase. The angular velocity experiences variation as long as the 
radiative core evolves and decouples from the convective enve¬ 
lope. In this process, all stars reach a common evolution after 
~ 1 Gyr. In our approach, we therefore propose to separate a syn¬ 
thetic parent sample into groups I and II (see Sect. O, with ages 
younger and older than 1 Gyr, respectively. If the theoretical (un¬ 
biased) distribution for group I are known, then observational 
biases can be implemented artificially to produce our synthetic 
sample. However, such an unbiased distribution is difficult to de¬ 
termine. We therefore discuss possible models that can be super¬ 
imposed with group II. 

Figure |4] presents a histogram with the empirical distribu¬ 
tion of rotation periods for a compilation of PM S Sun-like ro¬ 
tating stars obtained by iGallet & Bouvi^ (l2013h . The stars be¬ 
long to 13 star-forming regions or open clusters with ages be¬ 
tween 1 Myr and 1 Gyr and were selected within a mass range 
of 0.9 to 1.1 Mq. The figure also shows two fits with simple 
models: a Gaussian function and an exponential function. These 
fits were defined fro m two basic assumptions. First, the empir¬ 
ical distribution from IGallet & Bouvied (1201 3l) can be basically 
interpreted as having two populations: one of fast rotators, with 
periods shorter than ~ 1 d, and another with a spread of rotation 
periods peaking at approximately 4 days and ranging from 1 to 
23 days. Second, the interest for a comparison with the solar ro¬ 
tation period lies in the region of decreasing probability (to the 
right of the peak), which is the part that can be superimposed 
with group II. 

Therefore, the Gaussian fit was obtained within the domain 
where Pmt > 1 d, that is, excluding the fast rotators. The fit peaks 



Fig. 4. Rot ation period distribution of the PMS Sun-like cluster stars 
selected by IGallet & BouvIe3 ( 120131) . These stars belong to 13 star¬ 
forming regions or open clusters with ages ranging from 1 Myr to 1 Gyr 
and masses between 0.9 and 1.1 Mq. The red dashed curve is a Gaussian 
fit with this distribution for periods greater than one day (see text). The 
blue dotted curve is an exponential fit with this distribution, where the 
origin of the exponential function was placed in the distribution peak 
(see text). 


at approximately 3.5 d and has a standard deviation cr ^ 2.6 d. 
A group of fast rotators was not identified in the CoRoT and Ke¬ 
pler Sun-like samples (see Sect. IS; thus, for our applications, 
no cutoff was needed at short periods to compute Gaussian fits. 
This most likely occurred because the period range of fast ro¬ 
tators coincides with many other phenomena, such as pulsation 
and binarity; thus, more data than simply photometry data are 
needed to distinguish variability signals within this range. The 
exponential fit was designated as an approximation for the re¬ 
gion of decreasing probability, in this case, for > 4 d. This 
fit, defined with an origin at 4 d, has an exponential decay rate 
of approximately -0.4 d '. _ 

Because the sample of iGallet & BouvieJ (1201 3h is a compi¬ 
lation of different works, a complex mixture of biases certainly 
affected this distribution. As a consequence, the precise unbiased 
distribution of group I is difficult to determine empirically. De¬ 
spite these facts, we represent the distribution of group I by the 
two functions considered in the fits to test their superpositions 
with group II. For each sample analyzed in the results (Sect. |5]l, 
different Gaussian and exponential functions were tested with 
free parameters to fit the corresponding data, which helped us to 
interpret the full P^ot distributions. 

5. Results and discussion 

In this section, we present the main results obtained from the 
subsamples considered in this work. For the CoRoT sample, we 
first present (Sect. 15.11) an evolutionary study of the full DMS 
sample before analyzing the Pmt distribution of the DMSq sam¬ 
ple. This study is useful for validating the period measurements 
of DM13 as being rotational and allows for an initial interpre¬ 
tation of the Prot empirical distribution in the HR diagram. In 
particular, there is a clear deviation between the observed Prot 
average of Sun-like stars and the solar value that is caused by 
biases (see below). Explaining this effect is the main motivation 
of this work. 

We then analyze the empirical Prot distribution of each Sun¬ 
like subsample and consider synthetic biased distributions to un- 
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Table 1. Best-fit solutions. 


Parameter “ 

DMSo 

Sample name 
MHo* 

RHo 

Pifd) 

5.5 +0.6 

12.3 ±2.0 

11.3 ±0.8 

0-, (d) 

2.7 ±1.4 

6.4 ±2.3 

5.4 ±1.3 

PI (%) 

84 ±12 

64+12 

72+14 

Pii (%) 

16 ±12 

36 ±12 

28 ±14 

Det. par. 

SIN 

w 

A (mmag) 

thres 

1.00 ±0.65 

0.25 ±0.08 

1.78 ±0.61 

Dks 

0.058 

0.034 

0.025 

Pks (%) 

5E 

0.65 

20 


Notes. 

“ P] and CT] refer to the normal distribution peak and standard devia¬ 
tion for group I, Pi and pu depict the proportion of stars in groups I 
and II with respect to the entire Sun-like sample, Del. par. and thres 
give the detectability parameter considered in the simulation and the 
corresponding threshold value obtained in the fit, and £)ks and Pks are 
the resulting K-S distance and K-S probability. 

For the MHo sample, the detectability parameter w was simulated 
from a calibration with the amplitude A for an approximate solution. 

Pks increases with decreasing Dks and decreases with increasing the 
sample sizes. Hence, the relatively low number of DMSo-sample targets 
produces a high probability even though the MHo and RHq samples 
have lower Dks values. 


derstand the observational biases. Because the synthetic samples 
were produced using models that fit the solar evolution without 
bias well, the solar rotation normality can be discussed from a 
comparison between the actual and modeled distributions. For 
an overview of these results. Table [1] summarizes the best-fit so¬ 
lutions obtained for the three Sun-like samples considered in this 
work. Details are provided in Sects. 1572145.41 

5.1. Evolution across the HR diagram and Aot validation for 
the CoRoT sample 

A first important result obtained from the DMS sample is the 
period distribution across the HR diagram. Figur e |5] shows this 
distrib ution plotted with theoretical tracks from lEkstrom et al.l 
(l2012h . The figure also shows a diagram of periods as a function 
of Teff. This rich combination of photometric periods, spectral 
parameters, and theoretical predictions encompasses an ample 
set of cross-checking features for our physical interpretations. 

Figure |5] (bottom panel) strongly supp orts the self¬ 
consis tence between the models considered by lEkstrom et al.l 
(I 2 OI 2 I) for rotating stars and the DM13 period values. The to¬ 
tal independence of the models and these observations means 
that the phys i cal in gredients considered in the models of 
lEkstrom et al.l (12012h are compatible with observations. At the 
same time, with the new data used here, this result validates the 
DM13 periods (previously considered as rotational candidates) 
as being rotational. Although false positives may be present in 
this sample, they would clearly be a negligible fraction. 

Interestingly, the DM13 sample illustrates that this CoRoT 
subset is mainly distributed around the theoretical track repre¬ 
senting the evolving Sun. This distribution is a result of the selec¬ 
tion methods described in DM13, which provided a subsample 
characterized by well-defined semi-sinusoidal photometric vari¬ 
ations (see Sect. 0. This suggests that the semi-sinusoidal signa¬ 
ture defined in DM13 is a special characteristic of late-type low- 
mass stars, with a high occurrence for Sun-like stars at different 
evolutionary stages. This may have important implications in the 
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Fig. 5. Top panel: HR diagram for the CoRoT ro tating candidates o f 
DM13 combined with the physical parameters of ISarro et all ( I2013I) . 
namely, the DMS sample. Circle sizes and colors repre sent the period 
interva ls shown in the legend. Theoretical tracks from lEkstrom et ahl 
ll2012h are plotted as solid lines for different masses, identified by corre¬ 
sponding labels. The colors represent theoretical rotation periods within 
the same intervals as in the legend. The typical Teff and log g uncertain¬ 
ties are illustrated by the error bars. The Sun-like sample is denoted 
within the ellipsis around the solar location. Bottom panel: photomet¬ 
ric period as a function of effective temperature for the same sample 
as the top panel. Th e red solid lines represent theoretical tracks from 
lEkstrom et all ( 120121) for four different masses. The Sun is denoted in 
both panels by the solar symbol. 


study of chromospheric activity and spot behaviors in rotating 
stars and in their relation to the Sun. Regarding the temperature- 
period diagram, the Sun location is sparsely populated, and many 
more stars are observed with short periods around the solar tem¬ 
perature. Biases must be considered when determining whether 
this represents a trend for the Sun rotating more slowly than Sun¬ 
like stars. 
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Fig. 6. Rotation period distribution of a subset of CoRoT stars with 
and log g similar to those of the Sun (the DMSo sample). The origi- 
nal sample was obtained from DM13 combined with lSarro et alj l l2013h 
(the DMS sample). The black histogram shows the actual distribution. 
The blue and red histograms and the dotted curve correspond to a mod¬ 
eled distribution of the actual data where biases were simulated. Blue 
depicts a group of young stars (group I), red depicts a group of MS and 
evolved stars (group II), and the dotted curve is the sum of these two 
group distributions. The latter, which is the full distribution, is the syn¬ 
thetic data considered to fit with observations. For illustrative purposes, 
the levels of the modeled distributions were adjusted to make the total 
number of synthetic stars similar to that of the actual sample. 


5.2. Prot distribution of the CoRoT sampie around the 
iocation of the Sun in the HR diagram 

Figure|6]shows the period distribution for a selection of the DMS 
sample within an elliptical region around Teff = 5772 + 300 K 
and logg = 4.44 + 0.4 dex, as depicted in Fig |5] (top panel), 
namely, the DMS© sample. The distribution is asymmetric, with 
a clear peak around four days and a long tail toward long peri¬ 
ods, and it can be interpreted using a separation into two groups. 
We do not claim any evidence for a bimodal distribution, but 
the observations are compatible with it. Based on this assump¬ 
tion, there is a high occurrence of young stars, corresponding 
to group I, compared to group II. The reason why considerably 
fewer stars are observed with an increasing period is explained 
by a detecta bility bias. Thi s bias has its origin in the relations 
identified by iMessina et al.l (1200 ll l2003h . where the variability 
amplitude decreases as long as its period increases. We provide 
a quantitative analysis of the rotational modulation detectabil¬ 
ity (related with amplitude) to understand how actual samples of 
photometrically rotating stars can be biased. Thus, our method 
is an important contribution to analyze P^ot by co mparing actual 
biase d samples with theory. Other authors (e.g., lAigrain et al.l 
1201 5h have recently investigated methods to provide better con¬ 
straints on the rotational modulation detectability. 

For the synthetic sample, we applied a normal distribution to 
generate the rotation periods of group I, as discussed in Sect. H] 
For this sample, the synthetic distribution fit observations well 
by centering at approximately 5.5 + 0.6 d with a standard devi¬ 
ation of approximately 2.7 + 1.4 d (see Table [TJ. These values 
are compatible with those computed in Sect. |4] for a compiled 
sample of PMS stars. Exponential decays were also tested, but 
they provided similar fits when considering the full synthetic 
distribution. For group II, the rotation periods were obtaine d 
from interpolating theoretical grids from lEkstrom et aP (1201 2h . 
as explained in Sect. The proportions p/ and p// segregated 
by the 1 Gyr isochrone were 84% and 16%, respectively, which 



Fig. 7. Rotation period distribution of a newer subset of Kepler stars 
with Teff and logg similar to those of the Sun (the MHp sa mple). The 
origi nal sample was ob tained from iMcOuillan et alj ( |2014) combined 
with iHuber et alj(l2014h . The symbols are the same as in Fig.[^ 


fits the observed distribution without any additional adjustment. 
Thus, the full superposition of groups I and II in Fig. |6] corre¬ 
sponds with the observed distribution within the limitations of 
the CoRoT data. We cannot infer any conclusion about the solar 
rotation normality from this CoRoT sample because the data are 
too sparse around the solar period. However, the compatibility 
between the synthetic and observed distributions illustrates that 
the synthetic sample was generated based on reasonable assump¬ 
tions and could be a start to validate our approach. In addition, 
this result provides a basic explanation for the apparent abnor¬ 
mality of the solar rotation compared to Sun-like stars, as de¬ 
scribed in Sect. O The Sun-like samples obtained in this work 
cover very many young stars. This means that the many targets 
that rotate faster than the Sun in our Sun-like samples are most 
likely young objects. 

The CoRoT sample analyzed above was limited by a rela¬ 
tively low number of stars in comparison to the Kepler samples 
because of at least two reasons. First, based on Fig. ^ CoRoT 
covers a V^ag range of ~12-18 mag and noise levels of ~600- 
10,000 ppm, while for Kepler the Vmag range is of ~8-18 mag 
and the noise levels of ~30-5000 ppm. As such, Kepler has 
a higher photometric sensitivity and can observe brighter stars 
than CoRoT, which means that Kepler has more potential to de¬ 
tect micro-variability signals (such as the rotational modulation) 
for many stars. Second, the catalog of physical parameters ob¬ 
tained bv ISarro et al.l (l2013h (which we combined with the cata¬ 
log of DM13) has relatively few stars (6,832 sources), whereas 
the catalog of iHuber et al.l (12014h provides physical parameters 
for almost all Kepler targets (196,468 stars). Kepler therefore 
has considerably more data for Prot measurements and for phys¬ 
ical parameters and hence provides more reliable results, as pre¬ 
sented below in Sects. 15731 and l5 .41 


5.3. Testing a recent Kepler sampie 

As a test. Fig. Q shows the period distri bution of a Sun-like Ke- 
pler sample of rotating stars obtained bv lMcOuillan et 


combined with the physical parameters of IHuber et al. 
namely the MHq sample. This stellar set was selected within an 
elliptical region around the solar location in the HR-diagram, 
as illustrated in Fig. [1] (bottom panel). This is considerably more 
populated than the CoRoT sample and could be selected closer to 
the solar metallicity, specifically -0.2 dex < [Fe/H] < 0.2 dex. 
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Therefore, it may allow for a more refined analysis of the rotation 
period distribution based on the synthetic calculations. However, 
for this case, simulating the detectability parameter w properly 
was particularly challenging; this limits our interpretaion of the 
MH 0 sample. We consider this case for an overview of the results 
and to discuss the importance of simulating the detectability pa¬ 
rameter well. For the RHq sample (see Sect. l5.41 i. we were able 
to properly simulate a detectability cutoff, and we obtained an 
optimal fit between simulation and observations. 

Overall, the Prot distribution also has a bimodal shape and 
can be decomposed into groups I and II, as performed in the 
DMSq sample. In this test, the proportions pj and pn initially 
set up by separating the parent sample into the regions below 
and above the 1 Gyr isochrone (which provided p/ = 81% and 
pii - 19%) did not fit the observations well. A better fit was ob¬ 
tained by adjusting these proportions to 64% and 36%, respec¬ 
tively. This handle is a fair adjustment related with the period- 
amplitude distribution and with the line that separates the two 
groups (see Sect. [3]l. The simulation of group I was simply de¬ 
fined as a normal distribution, which agreed with observations by 
peaking at approximately 12.3 + 2.0 d with a standard deviation 
of approximately 6.4 + 2.3 d. Hence, this distribution is centered 
at longer periods and has a broader spread than in the CoRoT 
sample. This difference may be related to the distinct galactic re¬ 
gions and fields of view (FOV) of the CoRoT and Kepler stellar 
sets, as analyzed in Sect. EH Exponential decays also provided 
similar fits, which did not improve the results. 

Regarding the superposition of group I with group II, the 
decay of the synthetic distribution to the right of the second 
peak reaches from slightly steeper to smoother than in the ob¬ 
served distribution. As a result, the synthetic sample presents a 
small excess of stars for periods somewhat above the solar value 
in comparison with observations. The discrepancy might be re¬ 
duced if practical difficulties of measuring long periods were 
implemented in our approach. To check for this possibility, we 
tested convolving the synthetic distribution with a linear decay 
that starts in a 100% probability at P^ot - 0 and decreases with a 
coefficient that was set as a free parameter. This function is com¬ 
patible with the recovery fraction as a function of the number of 
observed variability cycles, as analyzed in Sect. 2.2.2 of DM13. 
In this test an excess remained at long periods, so this contribu¬ 
tion could not be sufficient to adjust the excess. Therefore, such 
an alternative ingredient would not produce a degeneracy in the 
fit solution and the discrepancy is likely due to our limitations in 
simulating w. This conclusion is discussed further in Sect. 15.41 
where the simulation of Ryar was easily approximated to the A 
numerical values. 

The discrepancy occurred in Fig. Q illustrates the general 
effect observed in the distribution if a proper detectabil¬ 
ity cutoff is not performed in the synthetic samples. This effect 
is furthermore noticeable if the proportions pj and pu are ad¬ 
justed as free parameters without performing a detectability cut¬ 
off (see Sect.O. Overall, pj and pu can balance the peak levels 
of each group, whereas a detectability cutoff produces a fine- 
tuning of the distribution slope at long periods. In particular, the 
detectability thresholds obtained here as free parameters were 
substantially plausible (see Table [TJ; SIN — 1.00 and w — 0.25 
perfec tly match the cutoffs used by DM13 and lMcOuillan et ^ 
(l2014l) . respectively, whereas A - 1.78 mmag is compati b le with 
the threshold of Ryai = 3%o defined by iReinhold et al.l (l2013h . 
Despite the discrepancy in Fig. |7] the K-S test provided a non¬ 
null probability, and the synthetic sample fit observations rea¬ 
sonably well. Thus, this fact indicates the normality of the solar 
Prot compared to Sun-like stars based on the MH 0 sample. 



Fig. 8. Rotation period distribution of a subset of Kepler stars with 
Teff and logg similar to those of the Sun (namely, the RH p sample). 
The original sample was obtained from iReinhold et ^ ( l2013ll combined 
with iHuber et all ll2014h . The symbols are the same as in Fig.[^ 


5.4. tn-depth study of another Kepler sample 

Figure [8] presents the rotati on period distr i bution of a Sun¬ 
like subset of sources of IReinhold et al.l (12013h combined 
with iHuber et’aH (l2014l) . namely, the RH 0 sample. This is a 
highly populated stellar set, similar to the MHo sample, also 
selected within the elliptical region depicted in Fig. [T] (bottom 
panel) and within -0.2 dex < [Fe/H] < 0.2 dex. In addition, the 
detectability parameter of the RH 0 sample, Ry^, can be easily 
simulated because its numerical value is highly similar to that 
of the variability amplitude A. Therefore, this sample allows for 
the most detailed and complete analysis of the Prot distribution 
in the present work. 

The Prot distribution of the RH 0 sample is noticeably dif¬ 
ferent than that of the MH 0 sample. One reason is certainly the 
different time spans of the LCs (the RH 0 sample uses only Q3 
data, the MHq sample uses Q3-Q14 data), which facilitates long- 
period detection in the MH 0 sample. Despite the difference, the 
distribution is also bimodal and can be decomposed into groups I 
and II, as performed for the samples analyzed above. In the syn¬ 
thetic sample, when groups I and II were initially decomposed 
by the regions below and above the 1 Gyr isochrone (giving pj 
and pii of 81% and 19%, respectively), the Prot distribution fit 
the observations reasonably well. However, the best fit was ob¬ 
tained by adjusting p/ and pu to 72% and 28%, respectively. 
As suggested in Sect. H) a normal distribution was defined for 
group I, which agreed with observations by peaking at approxi¬ 
mately 11.3 + 0.8 d with a standard deviation of approximately 
5.4 + 1.3 d. Hence, this distribution is compatible with the RHq 
sample. Exponential decays also provided acceptable fits and 
produced similar behaviors, which did not improve the results. 

Eor the entire P^t distribution of this sample, a noteworthy fit 
was obtained because the simulation was able to properly repro¬ 
duce the bias caused by the detectability parameter Pyar- Indeed, 
we verified that for the RH 0 sample, if S jN was used as the 
detectability parameter instead of A, a discrepancy would occur 
between the simulation and observations in the distribution de¬ 
cay around the solar period, similar to the problem of the MH 0 
sample. Therefore, an appropriate simulation of the detectability 
parameter certainly contributes to an optimal fit with observa- 
ti ons. Because the synth etic distribution has the theoretical grid 
of lEkstrom et al.l (12012l) as a basic ingredient, which unbiasedly 
fit the Sun and Sun-like stars, the good fit indicates that the solar 
Prot is normal with regard to Sun-like stars. 


Article number, page 10 ofll3l 









































I. C. Leao et al.: Rotation period distribution of CoRoT and Kepler Sun-like stars 


Compared to the other CoRoT and Kepler samples, the pe¬ 
riod distribution of group I is compatible with the RHq sample, 
and centered around a longer period and with a broader spread 
than in the DMSq sample. This supports the possible relation of 
the group I distribution with the stellar population, as mentioned 
in Sect. 15.31 Indeed, the CoRoT fields lie toward the center and 
anti-center of the Galaxy, covering very many stars located near 
the Galactic plane. This region is expected to be m ore populated 
by yo ung stars than at high Galactic latitudes (e.g.. lMiglio et al.l 
l2013h . The CoRoT fields are also small, covering together a re¬ 
gion of ~8 square degrees. In contrast, the Kepler field is consid¬ 
erably larger (~10 5 square degrees) and covers a broad r ange of 
Galactic latitudes dMiglio et alJ2013llGrrardi et alJ201^ . These 
facts suggest a higher proportion and a narrower distribution of 
young stars in the CoRoT sample than in the Kepler sample. To 
test this hypothesis, we used the TRILEGAL code to obtain a 
sampling of the CoRoT and Kepler fields, covering nearly the 
same regions of their actual coordinates. The limiting magni¬ 
tude was chosen to be V = 16 mag for the CoRoT fields and 
Kp - 17 mag for the Kepler field. Then we restricted the Tefr, 
logg, and {MjH} (assumed to be ^[Ee///]) ranges according to 
our own selection of Sun-like stars. Based on this test, the age 
distribution of the CoRot and Kepler samples agrees well with 
the Prot distributions presented here, where the CoRoT fields in¬ 
clude a greater proportion of young stars than older stars. In the 
Kepler field simulation, the simulation is much broader, with a 
considerably fewer young stars than in the CoRoT sample. 

Despite this difference, either the CoRoT or Kepler fields 
comprise a high proportion of young objects. Overall, 64-84% 
of the Sun-like samples analyzed here are younger than ~1 Gyr. 
This result (restricted to Sun-like stars) is re a sonab le if com¬ 
pared with a related estimate from iMatt et alJ (l2014l) . valid for 
the whole sample of Kepler sources, based on a relatively simple 
stellar rotation model. According to these authors, ~95% of the 
whole sample of Kepler field stars with Prot measurements are 
younger than ~4 Gyr. 

In general, our study shows that the CoRoT and Kepler data 
can be well interpreted in terms of two populations, one young 
(group I), one older (group II), and with a different ratio between 
them. The proportions pi and p// for the CoRoT and Kepler 
samples considered here can be well interpreted by simulating 
the experiment sensitivity and biases. At least on the basis of 
the present data, ther e is no clear reason to believe that the bi¬ 
modality reported by iMcOuillan et al.l (l2013ai l2014l) can be ex¬ 
plained by selection biases. As such, we support the suggestion 
of these authors, with some addition. The bimodality is possibly 
explained by a “natural bias” (i.e., not dependent on measure¬ 
ment sensitivity, but related to evolutionary duration) consider¬ 
ing that two stellar populations experience evolutionary phases 
that have very different durations. These populations are, in fact, 
the groups I and II defined here. Therefore, the detectability bias 
is not the reason of the bimodality, although it influences the 
proportions p/ and p//. 

Of course, the normal distribution adopted here for group I 
is an approximation that may have a number of implicit bi¬ 
ases. In contrast, the distribution of group II is based on a more 
elaborate set of theoretical model s and assumpti o ns, in cluding 
the theoretical grids provided by lEkstrom et al.l (l2012h . mod¬ 
eled par ent distributions, the empir ical period-amplitude rela¬ 
tions by iMessina et al.l (1200iL l^03h . the selection by SjN, A, 
or w (which mimic detectability biases of the rotational modula¬ 
tion), and final filtering by Sun-like characteristics, all described 
in Sect. [2 Hence, the set of ingredients we presented allows us 
to study solar rotation by comparing actual and biased observa¬ 


tions with theoretical models, within a precision that depends on 
the current data and model limitations. More refined studies can 
be developed in the future with this method. In particular, a po¬ 
tential improvement of the current results would be possible if 
spot-induced variability amplitudes were included in the stellar 
evolution codes with rotation. 

6. Conclusions 

We presented for the first time a quantified method with which 
the empirical distribution of Prot for Sun-like stars can be ex¬ 
plained by considering the detectability of rotational modulation, 
which is related to the variability amplitude. The interpretations 
were based on a comparison between synthetic and empirical 
subsamples of CoRoT and Kepler field stars with physical pa¬ 
rameters similar to those of the Sun. 

Eor this purpose, we combined the public data of rotation pe¬ 
riods with physical parameters for three stellar samples. The first 
sample was obtained from a CoRoT catalog with period mea¬ 
surements given by DM13, combined with spectroscopic phys - 
ical parameters log.g and determined bv ISarro et al.l (1201 3l) . 
namely, the DMS sample. The selection of Sun-like stars was 
obtained within an elliptical region around the solar location 
in the HR diagram. Eor CoRoT stars, the subsample was en¬ 
closed by Teff = 5772 + 300 K and logg = 4.44 + 0.4 dex, 
namely, the DMS© sample. The two other samples were se- 
lected from Kepler c atalo gs of rotation periods provided by 
iReinhold et al.l (1201 3h and iMcOuillan et ^ (l2014t) . each com¬ 
bined w ith an improved c ompilation of physical parameters 
given bv iHuber et alJ (l2014l) . The Sun-like subsamples were se¬ 
lected within -0.2 dex < {FejH} < 0.2 dex and within an ellipsis 
described by — 5772 + 170 K and logg = 4.44 + 0.2 dex, 
namely, the RH© and MHq samples. 

To explain the observational biases, synthetic samples were 
generated based on the parent samples of the CoRoT and Ke¬ 
pler fields. The synthetic stars were separated into two groups, 
I and II, which are those located below and above the 1 Gyr 
isochrone in the HR diagram, respectively. These represent 
a group of young stars and a group composed of MS and 
evolved stars, respectively. Eor group I, synthetic rotation pe¬ 
riods were randomly generated following two possible distribu¬ 
tions compatible with PMS Sun-like stars: normal or an expo- 
nenti al decay from the pe ak. Eor group II, theoretical periods 
from lEkstrom et akl (1201 2l) were set up according to the HR di¬ 
agram locati on of the synthetic stars. Ne xt, period-amplitude re¬ 
lations from iMessina et all (1200 ll [20031) were used to aggregate 
amplitude values to the synthetic samples, and random noise lev¬ 
els were added to these samples. Then, the main bias related to 
the detectability of the rotational modulation was applied to the 
synthetic stars by selecti ng them with a thresho ld in the variabil¬ 
ity SIN, A, or w (from IMcOuillan et al.ll2014l which was par¬ 
ticularly difficult to simulate) according to each observed sam¬ 
ple. Einally, Sun-like subsamples were obtained by selecting the 
targets located within the ellipses described around the solar lo¬ 
cation in the HR diagram from the synthetic samples. The el¬ 
lipses were defined within the same regions as for each one of 
the three observational subsamples of Sun-like stars considered 
here, namely, the DMS©, RHq, and MHq samples. 

The distribution in the HR diagram of the DMS sample has a 
large number of stars with low masses that are particularly sim¬ 
ilar to that of the Sun. This sample was obtained by identify¬ 
ing a well-defined semi-sinusoidal signature in the light curves 
without knowledge of their physical parameters, apart from pho¬ 
tometric magnitudes and colors. This sample indicated that the 
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semi-sinusoidal signature, defined by DM13, provides a good 
selection of low-mass spotted rotating stars, with a high occur¬ 
rence for stars with masses of ~1 Mq, which may have important 
implications in the relation between rotation and chromospheric 
activity for evolving Sun-like stars. 

The rotation period distribution of the DMS© sample is not 
clearly bimodal, but it hts our bimodal model, thus showing a 
high occurrence of young stars. The proportions of stars belong¬ 
ing to groups I and II with respect to the total number of ob¬ 
jects in the synthetic sample, namely, p/ and p//, are of 84% and 
16%, respectively. These proportions ht the observations when 
simply separated by the 1 Gyr isochrone without any comple¬ 
mentary adjustment. In addition, the full distribution of the syn¬ 
thetic sample, that is, the sum of groups I and II, exhibited a 
close ht with the observed distribution. This shows that separat¬ 
ing the two groups by using the 1 Gyr isochrone is reasonable. 
Because of the relatively low number of objects in the DMSq 
sample (175 targets), we cannot make conclusions regarding the 
solar rotation normality with respect to its Sun-like counterparts. 
Therefore, this result helps to validate the methods used to pro¬ 
duce our biased synthetic samples. Further in-depth analyses can 
be performed from the Kepler data. 

For the MHq sample, the period distribution is more seem¬ 
ingly bimodal and has considerably more stars (2,525 targets) 
than the DMS© sample. For this case, the proportions p/ and p// 
originally did not ht the observations well and was adjusted to 
hnal proportions of 64% and 36%, respectively. With this adjust¬ 
ment, the modeled distribution ht the observations well, although 
it exhibited a noticeable discrepancy in the distribution decay 
that traverses the solar rotation period. However, near the so¬ 
lar period, the synthetic sample hts the observations well, which 
suggests that the solar rotation period is normal with respect to 
Sun-like stars. 

The period distribution of the RHq sample also shows a bi¬ 
modality and has a high number of stars (1,836 targets), sim¬ 
ilar to the MHq sample. The proportions p/ and p// originally 
ht the observations well, although the best ht was obtained by 
adjusting these parameters to the hnal values of 72% and 28%, 
respectively. With this adjustment, the full period distribution of 
the synthetic sample provided an extremely close ht with obser¬ 
vations. This fact strongly supports that the solar Prot is normal 
with respect to Sun-like stars. 

Overall, the main conclusions of this work can be summa¬ 
rized as follows: 

i. The synthetic samples were able to explain the observed 
Prot distributions by resolving them into two different pop¬ 
ulations, group I for young stars and group II for main- 
sequence and evolved stars, which generally produce a bi¬ 
modal arrangement. 

ii. The central peak and the width of the group I P^ot distri¬ 
bution may be related with the stellar population within a 
certain FOV. 

iii. Reasonable to optimal hts of the synthetic samples with the 
observed data were obtained especially by adjusting the bi¬ 
modal peak levels with p/ and p// and by rehning the distri¬ 
bution slope at long periods with an appropriate detectabil¬ 
ity threshold for each Sun-like sample. 

iv. Several tests indicated that the parameters considered in 
these hts (see Table [T]) provide a simple description of the 
Prot distributions without degeneracy. 

V. A detectability threshold is mandatory to properly ht syn¬ 
thetic with observed data for the CoRoT sample, while it 
improves the hts for the Kepler samples, but with a weaker 


effect. This suggests the Kepler samples are less biased by 
the variability amplitude than the CoRoT sample. 

vi. Best hts suggest a high number (64-84%) of young objects 
in the Sun-like samples, which explains the short Prot values 
of Sun-like stars on average compared to the solar value. 

vii. The global agreement between the synthetic and observed 
samples suggests the normality of the solar Prot, at least 
within the current data accuracy and model limitations. 

Therefore, the method presented in this work allowed us to con¬ 
strain observations on theory by considering biased samples with 
Prot measurements and was particularly useful to analyze held 
Sun-like stars. In addition, the following perspectives can im¬ 
prove our results: 

i. A substantial improvement of the synthetic samples will be 
possible once stellar evolution models can predict the photo¬ 
metric amplitude variations produced by rotating star spots. 
Such a theoretical prediction needs a description of the mag¬ 
netic activity at the surface of stars, and the interpretation of 
the observed relation between P^t and A is one of the ob¬ 
servational facts that can help in this subject. 

ii. Observed data with more accurate physical parameters will 
provide a more rehned selection of Sun-like stars and a more 
reliable separation of groups I and II. Accordingly, future 
spectroscopic observations, s uch as those to be collected 
with the Gaia-ESO SurvejQ (iGilmore et al.ll2()T^ . will al¬ 
low a better verihcation of our results. 
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